Passive Sliders on Growing Surfaces and (anti-) Advection in Burger's Flows 
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We study the fluctuations of particles sliding on a stochastically growing surface. This problem can 
be mapped to motion of passive scalars in a randomly stirred Burger's flow. Renormalization group 
studies, simulations, and scaling arguments in one dimension, suggest a rich set of phenomena: If 
particles slide with the avalanche of growth sites (advection with the fluid) , they tend to cluster and 
follow the surface dynamics. However, for particles sliding against the avalanche (anti-advection) , 
we find slower diffusion dynamics, and density fluctuations with no simple relation to the underlying 
fluid, possibly with continuously varying exponents. 



PACS numbers: 68.35.Rh, 02.50.Ey, 64.60.Ht 

The advection of a passive scalar, such as dye, a tem- 
perature field, or dust particles, in a turbulent flow is a 
scientific problem of great interest, since the fluctuations 
of the scalar field appear to have no simple relation to the 
statistics of the flow [Q . This holds not only for fully de- 
velopped turbulence, as emerging from the Navier-Stokes 
equation, but even for the much simpler case of a random 
Gaussian velocity field which has power-law correlations 
in space and is delta-correlated in time ||. 

In this article, we introduce and study a related prob- 
lem of particles passively sliding on a fluctuating land- 
scape, such as a growing surface. To exploit its mathe- 
matical similarity to the better known problem of advec- 
tion of a passive scalar, we start by considering a velocity 
field v(x, t) which follows the Burger's equation, given by 

dv 

— + X(v ■ V)v — vV 2 v + V(h{x, t). (1) 
at 

We also impose the additional condition 



-V/i, 



which ensures a gradient flow with no vortices. This con- 
dition is preserved by the random stirring with the gra- 
dient of a Gaussian white noise with correlations 

(( h (x, t)C h @', f)) = 2D h 6 d (x - x')5(t - t'). 

The noise term also prevents the formation of shock 
waves. Shock waves occur in the absence of a stirring 
force, or with stirring correlated over large spatial scales 
f|. In contrast to the Navier-Stokes equation, Burger's 
equation has no pressure term, and the fluid is there- 
fore compressible. Our use of Burger's equation is not so 
much motivated by describing a particular physical situa- 
tion, but because its flows are simpler than Navier-Stokes 
turbulence, which is a yet unsolved problem. When ex- 
pressed in terms of the scalar variable h(x,t), Eq. (Q) 
becomes the Kardar-Parisi-Zhang (KPZ) equation [Q for 



interface growth, 
dh 



nt = vV 2 h+-(Vh) 2 +( h (x,t), 



(2) 



where h(x, t) is the height of interface at position x at 
time t. Fluctuations of the surface height obey the scaling 
relation || 



([h(x,t) - h(x' ',t')f) = \x-x" | 2x . 



t-t' 



\x — x 



(3) 



In one dimension, the only case that we shall consider 
here, the roughness exponent is x — an( i ^ ne dy~ 
namic exponent is z = 3/2. These values are exact due 
to a fluctuation-dissipation condition, and the Galilei in- 
variance of Eq. (Q) (|] . 

A particle moved by the flow is characterized by the 
Langevin-equation 



dx 
~dt 



av - 



•£»(*), 



(4) 



where the white noise ( x describes the effects of the tem- 
perature and has a correlation function 

(Cx(i)C*(*')> = 2K<5(t-i'). 

For a — A, the particle is advected in the same way as 
an element of the fluid. For < a < A, the particle is 
slower than the neighboring fluid elements, while a nega- 
tive a describes a particle moving against the flow (anti- 
advection). In the language of the fluctuating interface 
(with v = — V/i), Eq. (||) describes a particle sitting on 
the surface and sliding with a velocity proportional to 
the local slope, for instance under the influence of grav- 
itation. (This equation also describes the dynamics of a 
single domain wall separating two phases 0.) The cor- 
responding equation for the particle density p(x, t) is 



dp 
dt 



nV 2 p + aV(pVh) + C P (x,t), 



(5) 



2 



where we have again added a white noise, with 

(( p (x, t)( p (x ', 0) = -2D p V 2 6 d (x - x')S(t - t 1 ) . 

In this paper, we consider only nonintcracting particles 
that move independently of each other. The case of many 
particles with hard-core repulsion is treated in Rcf. 

The difficulty in analysis of the stochastic Eqs. (|2| 
arises from the non- linear couplings. A perturbative 
treatment of the non-linearities results in divergences 
that have to be controlled, e.g. by a renormalization 
group (RG) treatment. Such a treatment is well-known 
for the KPZ equation (||), and can be extended to the 
advection equation (^|). In particular, at the one loop 
order, the reduced parameters a = a/X and k — k/v, 
satisfy the RG flows 



da 
~dl 
dk 
~dl 
dD p 
dl 



2{a-a 2 ) 



k 
D 



(1 + fc) 2 
a 2 (3-fc) 



k(l + k) 2 



k(l + k) 



(6) 



Similar to Eq. (||), the exponents \ p and z p describe the 
scaling of particle density fluctuations via 



([p(x,t)-p(x',t')] 2 )~\ x - x '\ 



2Xp 



t-t' 



(7) 

Standard dynamic RG calls for a single exponent z for 
the relative scaling of time and space. The reason for 
including a distinct z p for the dynamics of particles shall 
be explained below. Other than this difference, the above 
RG equations are a special case of the flows derived in 
Ref. Q for a polymer drifting through a random medium. 

Equations (^) show quite distinct behavior depending 
on the sign of a — a/X. While a > is the natural 
choice for advection by a fluid, both signs are possible 
for particles sliding on a fluctuating surface. Indeed, we 
first encountered a < (anti-advection) in connection 
with domain walls in a model of growing surfaces 0, 
inspiring the further simulations described below. For 
a > 0, the RG equations have an attractive fixed point 
at finite parameter values, with the standard dynamic 
exponent z — z p = 3/2. Using z p — 3/2, it is easy 
to show that the first two RG equations (||) yield the 
fixed point values a* = 1 and the k* = 1, while the 
third sets \p = 3/4. Note that for a = k = 1, the den- 
sity fluctuations share the Galilei-invariance of the fluid 
flows, and since symmetries are preserved under renor- 
malization, this fixed point could have been predicted 
without any calculations. However, this symmetry does 
not fix the scaling dimension x P ( m contrast to the val- 
ues of x and z), and we shall in fact show below that the 
observed value of Xp = 1/2 is different from the above 
one-loop prediction. The case a/X > is the one that 



is usually considered in related work, as for instance in 
Ref. |@[, which contains a paragraph on passive scalar 
advection in randomly stirred flows; in Ref. jl(| which 
discusses (among others) a situation that is equivalent 
to the advection of a single particle by Burger's flow; in 
Ref. [[ll] which focusscs on the world lines of particles, 
and in Ref. @. 

For a/X < and z p — z, the RG flows drive the pa- 
rameters a and k towards zero, such that to leading or- 
der a 2 /k = 1/6 (and \ p — 5/12). The vanishing of the 
parameter a under the coarse-graining action of RG sug- 
gests that on larger scales the particles are less able to 
perform their sliding motion. A potential interpretation 
of the lack of a standard fixed point is that the particles 
do not follow the fluctuating surface, and are instead gov- 
erned by a slower timescale. (A similar vanishing of pa- 
rameters under RG occurs for decoupled equations with 
distinct time scales.) Allowing for a distinct value of z p 
in the physically meaningful interval (3/2,2), does indeed 
lead to nonzero fixed points for a and k. Interestingly, the 
resulting exponent is non-universal, varying continuously 
along the fixed line a*(k). This conclusion appears to 
be supported by the simulations described below. How- 
ever, we shall make no quantitative claims based on this 
RG study, both because of the unconventional use of two 
timescales, and the "strong-coupling" nature of the fixed 
point, which precludes controlled perturbation expansion 
in any dimension. 

Next, we present the results of computer simulations. 
The growing interface is obtained from a restricted solid- 
on solid model, by adding "bricks" to a growing "wall" . 
Starting from a flat surface and imposing periodic bound- 
ary conditions, bricks are added to randomly chosen sites 
such that no overhangs are formed, and with the center 
of each brick atop the intersection of two bricks in the 
layer below. Such a growth rule leads to surfaces that 
are in the same universality class as the KPZ equation. 
The parameter A is negative in this model as the average 
growth velocity is less for a sloped surface. The advected 
particles are intially placed at random surface positions. 
Whenever a brick is added to the site of a particle, it is 
moved to the lower one of its two nearest-neighbor sites 
with probability (1 + K)/2, and to the higher nearest- 
neighbor site with probability (1 — K)/2. The parameter 
K is adjustable, but if not explicitly mentioned in the fol- 
lowing, it is set to K = 1. When both nearest-neighbor 
sites are at the same height, the particle is moved to the 
right or left with equal probability. The downhill charac- 
ter of these moves corresponds to a > in Eq. (||). The 
case a < is implemented similarly, moving particles 
uphill instead of downhill. 

After a time period long enough for establishing a sta- 
tionary state, the diffusion of the particles and the par- 
ticle densities were quantified. Figure [l] shows the mean 
square distance (x 2 (t)) ~ t 2 l Zp travelled by a particle as 
function of time, for different values of the simulation 
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FIG. 1: The mean square distance travelled by the sliding 
particles as function of time. The top curve is for a/A > 
(advection), while the others are for a/A < (anti-advection) 
with K = 1, 0.25, 0.125, 0.025 (from top to bottom). The 
dashed straight lines are power-law fits to exponents 2/z p , 
with z p equal to 1.5, 1.74, 1.83, 1.91, 1.98 from top to bottom. 



parameters. As suggested by the RG analysis, the parti- 
cles can follow the surface fluctuations for aj A > 0, and 
their diffusion is characterized by the same dynamical ex- 
ponent z p — 3/2 as the surface. For a/A < (particles 
moving downhill) , the diffusion exponent z p is larger, and 
it seems to be nonuniversal, i.e., it depends on the value 
of the parameter K. 

This qualitative difference between the cases a/A > 
and a I A < is also visible in the correlations of particle 
density. To examine this quantity, we divided the system 
into L/l bins of size / and counted the number of particles 
in each bin. Let us choose a particle at random, and 
denote by N(l, L) the number of particles that are in the 
same bin. This quantity is plotted in Fig. |^ as function 
of I, for a/A > 0. Using the scaling form of Eq. (]?]), we 
have 

N(l,L) ~ I^J ' dxp{x,t)p{0,t)dx^ ~ /2 Xp -i L 2-2 XPj 

(8) 

where the dependence on L follows from N(L, L) being 
the total number of particles in the system, which was 
chosen proportional to the system size L. The curves 
in Fig. H are in fact proportional to L, and become con- 
stant for larger I, both results consistent with \ p — 1/2 in 
Eq. (JsJ) . This value is different from the 1-loop RG result 
obtained above. The corresponding results for a/A < 
are shown in Fig. |^ (for K = 1). Since particles cannot 
fully follow the surface fluctuations, the proportion of 
particles within a given interval around a randomly cho- 
sen particle increases with increasing system size, and we 
measured a nontrivial scaling exponent Xp — 0.85. For 
smaller values of K, the scaling behavior is not so clear. 
The asymptotic scaling of the factor that maps curves 
for different L onto each other could not be determined 
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FIG. 2: Mean number of particles in the same bin of size I as 
a randomly selected particle, for a/A > and different L. 
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FIG. 3: Mean number of particles in the same bin of size I 
as a randomly selected particle, for a/A < and different L. 
The straight line in the bottom figure is a power law with 
exponent 0.7. 

reliably. It appears, however, that Xp is nonuniversal, be- 
coming larger with decreasing K (and increasing particle 
diffusion exponent z p ). 

Particle clustering, which is the most prominent fea- 
ture of our system, does not occur in Navier-Stokes tur- 
bulence, where the fluid is incompressible. Instead, fluid 
elements are streched and folded, leading ultimately to a 
homogeneous particle distribution, unless new particles 
are continuously added. 

In order to better understand the nonuniversal and 
nontrivial scaling for a/ A < 0, we measured directly the 
velocity of particles as a function of an external slope. A 
particular slope was imposed on the system by choosing 
periodic boundary conditions with a fixed height differ- 
ence between the two ends. The average speed of parti- 
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Finally, we shall try to compute the exponent y by 
relating characteristic time scales of surface and particle 
fluctuations. Consider a portion of the surface of size I, 
with a typical slope of order |V/i| ~ l~ x / 2 . For this slope 
to change by surface fluctuations, we should wait a time 
of the order t ~ l z ~ Z 3 / 2 . During this time the advecting 
particle will move against the slope by an amount 

leading to the scaling relation 



FIG. 4: The scaling function for the particle velocity as a 
function of slope and system size, for a/X < and K = 1. 
The straight line has the asymptotic slope 2y, with y = 0.14. 



cles was then evaluated as function of slope and system 
size. As expected, for a/X > the mean velocity is sim- 
ply proportional to the slope. For a/X < 0, however, the 
velocity shows a nontrivial scaling behavior of the form 

v = L- y VhC(Vh^/L), (9) 

with y ~ 0.14, as shown in Figure ^. The scaling func- 
tion C(x) is constant for small x and increases as x 2y for 
large x, such that v oc (Vh) 1+2y is independent of L. 
(Note that the data for L = 1024 show large fluctuations 
and seem sensitive to the choice of the random number 
generator, and therefore not all data points collapse with 
those for smaller L.) Just as for the density correlations, 
the data for smaller values of K did not show nice scaling 
behavior for the system sizes in our simulations. Equa- 
tion (||) implies that the factor a relating velocity to slope 
[see Eq. (0)] decreases as a function of system size, consis- 
tent with the RG prediction. However, while the 1-loop 
result suggests an inverse logarithmic dependence of a on 
the length scale (implying y = 0), the value of y found in 
the simulations is nonzero. 

The form of Eq. @ can be justified by the following 
argument: The growth rules for the surface imply that 
sequences of brick addition usually move from local min- 
ima in the uphill direction. The reason is that the addi- 
tion of a brick generates a potential growth site (where a 
brick can be added without generating overhangs) at the 
nearest uphill position. The particles that try to slide 
downhill are therefore faced with an upward avalanche 
of growth mounds of different sizes that hamper their 
downhill motion. (This is analogous to a moving fluid 
density fluctuation faced by a particle advecting against 
the flow.) The horizontal size of the largest mounds, I, 
scales as LC(Vh^/L), with C(x) being constant for small 



2(1 - y) 



(10) 



x and C(x) 



for large x. If we assume that the 



particle velocity is slowed down by some power of the 
size of the largest mounds, we obtain Eq. (|). The expo- 
nent y may then have a non-universal dependence on the 
manner in which the particle goes over these barriers. 



For y ~ 0.14, this gives z p ~ 1.74, which is the value 
measured in simulations of particle diffusion for K = 1. 

In conclusion, we have found surprisingly intricate scal- 
ing phenomena for passive sliders on stochastically grow- 
ing surfaces, and the equivalent problem of passive scalar 
motion in Burger's turbulence. In particular, there is 
a qualitative difference between the cases a/X > and 
a/X < 0, describing respectively passive advection with 
or against the flow. The former leads to clustering, while 
the latter yields non-trivial exponents. This conclusion is 
bolstered by RG analysis, simulations, and scaling argu- 
ments. While the simulations provide convincing indica- 
tion of a nonuniversal value of the exponent z p character- 
izing particle diffusion(for a/X < 0), the evidence is less 
compelling for as the density correlations. The absence 
of nice scaling for most values of simulation parameters 
may indicate slow crossovers, but then also leaves open 
the possibility that the exponent z p is universal even in 
the case a/X < 0. 
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